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Abstract. We consider a model mixture of hard colloidal spheres and non-adsorbing 
polymer chains in a theta solvent. The polymer component is modelled as a 
polydisperse mixture of effective spheres, mutually noninteracting but excluded from 
the colloids, with radii that are free to adjust to allow for colloid-induced compression. 
We investigate the bulk fluid demixing behaviour of this model system using a 
geometry-based density-functional theory that includes the polymer size polydispersity 
and configurational free energy, obtained from the exact radius-of-gyration distribution 
for an ideal (random-walk) chain. Free energies are computed by minimizing the 
free energy functional with respect to the polymer size distribution. With increasing 
colloid concentration and polymer-to-colloid size ratio, colloidal confinement is found 
to increasingly compress the polymers. Correspondingly, the demixing fluid binodal 
shifts, compared to the incompressible-polymer binodal, to higher polymer densities 
on the colloid-rich branch, stabilizing the mixed phase. 
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1. Introduction 

The physical properties of soft matter systems are often influenced by size polydispersity 
on mesoscopic length scales In colloidal suspensions, for example, variations 
in particle size can profoundly affect fluid-solid transitions 0. In polymer solutions, 
chain length (degree of polymerization) often has a broad distribution. In practice, 
distributions in colloidal particle radius and polymer chain length often can be narrowed 
by physical selection or special synthesis methods. However, even monodisperse 
polymer chains have a radius of gyration that fluctuates in response to the surrounding 
environment. In this sense, polymers in solution are fundamentally polydisperse. 

Polydispersity is especially relevant in mixtures of colloids and nonadsorbing (free) 
polymers, in which the polymer size sets the range of effective depletion-induced 
interactions between colloidal particles. A popular basis for modelling colloid-polymer 
(CP) mixtures is the Asakura-Oosawa (AO) model |3]. The AO model neglects 
fluctuations in chain conformations by approximating the polymers as effective spheres 
of fixed radius that are excluded from the hard-sphere (HS) colloids. The model 
further ignores polymer-polymer interactions by treating the polymer spheres as freely 
interpenetrating, an approximation at best justified in a theta solvent. 

Despite its simplicity, the AO model qualitatively captures thermodynamic phase 
behaviour observed in real CP mixtures - in particular, bulk fluid demixing into colloid- 
rich and -poor phases |S1 13 El • To address quantitative discrepancies with experiment, 
several recent studies have gone further to include polymer-polymer interactions - in 
both good solvents UHl HI] and poor [12] - and polydispersity in polymer chain 
length [13 CH US]. In all of these studies, however, the polymers are modelled as 
effective spheres of fixed size (or fixed size distribution). 

An issue that has not yet been widely explored is the influence of colloidal 
confinement on polymer conformations and implications for phase behaviour of CP 
mixtures. Polymers near surfaces and polymer-mediated interactions between surfaces 
have been examined by a variety of theoretical approaches, including scaling theory ^Hl 
[TTIIIH], integral equation theory jUJi^D], classical density- functional (DF) theory pT| . 
and field-theoretic methods |22]- A far more challenging problem is the influence on 
polymers of confinement to a fluctuating random porous medium created by colloidal 
particles in suspension. 

Numerical studies of colloids mixed with segmented polymer chains give some 
insight into polymer conformations in CP mixtures. Meijer and Frenkel |^ performed 
Monte Carlo simulations of hard spheres mixed with ideal (noninteracting) lattice- 
polymers. Dickman and Yethiraj |21] simulated mixtures of hard spheres and off- 
lattice polymer chains, modelled as freely-jointed "pearl necklaces" of hard spheres. 
Simulations such as these have revealed the potential importance of fluctuations in 
polymer size and shape in CP mixtures. 

Here we consider an intermediate model in which the polymers maintain 
spherical shape but can adjust their radius-of-gyration distribution to the bulk colloid 
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concentration. The constraint of fixed polymer size is tlius relaxed by endowing the 
effective polymer spheres with an internal degree of freedom that can respond to colloidal 
confinement. We study the phase behaviour of a model mixture of HS colloids and 
compressible polymers by means of a geometry-based DF theory. Confinement- induced 
adjustments in polymer size are permitted by including in the free energy functional the 
free energy of ideal polymers, derived from the exact radius-of-gyration distribution of 
a random- walk chain. 

The theory, when applied to bulk fluid phases, predicts both the demixing binodal 
and the polymer size distribution as a function of colloid and polymer concentrations. 
For sufficiently large polymer-to-colloid size ratio, the model displays bulk fluid demixing 
into colloid-rich and -poor phases, qualitatively similar to the behaviour of the AO 
model. Upon demixing, however, the relatively unconfined polymers in the colloid-poor 
(vapour) phase coexist with compressed polymers in the colloid-rich (liquid) phase. We 
find that the polymer radius of gyration can be significantly reduced, by up to 20 percent 
or more, in the colloidal liquid phase. Correspondingly, the fluid demixing binodal shifts 
relative to the incompressible-polymer binodal. 

Next, in Sec. we explicitly define the model system. In Sec. El we construct an 
appropriate classical DF theory that incorporates the conformational free energy of the 
polymer chains. Results for demixing phase diagrams and polymer size distributions 
are presented in Sec. IH Finally, we conclude in Sec. El 

2. Model 

The model we consider is a generalization of the AO model to the case of compressible 
polymers. Explicitly, we consider a mixture of hard colloidal spheres (species C), 
of monodisperse radius Rc, and non-adsorbing, linear polymer chains (species P), 
monodisperse in length, but polydisperse in radius of gyration, Rp, in a volume V 
at a temperature T near the theta temperature. The restriction to a theta solvent 
allows the polymer chains to be reasonably approximated as noninteracting random 
walks. The colloids interact via a hard-sphere pair potential: vcc{r) = cxo, if r < 2Rc, 
zero otherwise, where r is the centre-centre interparticle distance. When interacting 
with colloids, each polymer is assumed to behave as an effective hard sphere of radius 
equal to its radius of gyration: vcp{r) = oo, if r < Rc + Rp, zero otherwise. Finally, 
interactions between polymers vanish for all distances: vpp{r) = 0. 

To fix the polymer chemical potential, it is convenient to imagine a reservoir of 
pure polymer solution with which the system can freely exchange polymer (but not 
colloid) through a semi-permeable membrane. The reservoir also serves as a colloid-free 
reference state in which the polymer assumes an ideal radius-of-gyration distribution. 
Bulk fluid states are specifled by the mean colloid and polymer (reservoir) number 
densities, pc and pp (pp), respectively. As bulk thermodynamic parameters, we use the 
colloid packing fraction, rjc = (47r/3)p(7_R^, and effective polymer packing fractions, 
rjp = {4:7[/3)pp{R^g)^ and rjp = {An / 3) pp{RgY , where i?^ denotes the root-mean- 
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square (rms) radius of gyration in the reservoir, a quantity that is directly accessible 
in experiments. The reservoir polymer-to-colloid size ratio, Rg/Rc, provides a useful 
control parameter for tuning interparticle interactions, and thus thermodynamics. 

3. Theory 

3.1. Density Functional Theory 

To investigate thermodynamic properties of the model system, we focus on the 
Helmholtz free energy as a functional of the inhomogeneous density profiles: the colloid 
density, pci''^), and a continuum of polymer densities, {pp(r; Rp)}, indexed by Rp and 
normalized to the mean polymer density via J d^r dRp pp{r] Rp) = pp. It is 
convenient to separate the total free energy functional into three terms: 

F[pc(r), {pp(r; Rp)}] = [pc(r), {pp(r; Rp)}] + F,, [pc(r), {pp(r; Rp)}] 

/roo 
d\ dRp pp{v- i?p)/ehain(r; Rp', [pc(r)]), (1) 

where is the ideal free energy, is the excess free energy in the AO model 
generalized to polydisperse (but incompressible) polymer, and /cham(r; -Rp; [pc(r)]) is 
the local free energy of a compressible polymer chain. The third term in Eq. stems 
from the internal degrees of freedom of the chains and is formally equivalent to the 
contribution of an external potential acting on the polymers. The ideal free energy, 
associated with (centre-of-mass) translational and mixing entropy of the colloids and 
polymers, is given exactly by 

= j dV pc(r) [in (pc(r)A^) - I 

+ / dV dRp pp(r; Rp) [in (pp(r; Rp)kl) - l] , (2) 

where = l/k-oT., k-Q is Boltzmann's constant, and A,^ and Ap are the respective colloid 
and polymer thermal wavelengths. 

The excess free energy, arising from interactions, must be approximated. For this 
purpose, we adopt a geometry-based DF approach, which is immediately apphcable to 
multi-component systems. As a basis, we start from a recently-proposed DF theory [2Z] 
for a binary CP mixture in the AO model, the homogeneous limit of which is equivalent 
to the free volume theory of Lekkerkerker et al [Zj . Here we generalize this theory [2Zj 
to mixtures of monodisperse colloids and polydisperse polymers. Following previous 
work [23 1211 12H1 12n] , the excess free energy is expressed in the form 

/3Fe. = j d^x* ({n^(x)}, K(x)}) , (3) 

where the excess free energy density, $, is a function of a set of weighted densities for 
colloid and polymer species. The weighted densities are defined as convolutions, with 
respect to geometric weight functions of the actual density profiles and, in the case of 
polymers, integration over radius of gyration: 

n^; (x) = / dV pc{v) w,{^ - r; Rc), (4) 
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d\ dRp pp(r; Rp) w,{x - r; Rp). (5) 

We use standard fundamental-measure weight functions [23 I^E], u]^{r;R), 
V =0,l,2,3,vl,v2,m2, for spheres of radius R, given by 

w^iy- R) = 6{R - r), w,{r- R) = Q{R - r), (6) 



w,2(r; R) = 6iR - r) ^, ^^^{v- i?) = - r) ( ^ - ^ ) , (7) 



where r = |r|, 5{r) is the Dirac distribution, B(r) is the step function, and 1 is 
the identity matrix. Further, hnearly dependent weight functions are iy^(r;i?) = 
w^iv; R)/{A7rR), Wf){r; R) = w^{r; R)/R, and w^i(r; R) = w^2(r; i?)/(47ri?). The weight 
functions have dimensions (length)''"^ and differ in tensorial rank: Wf),Wi,W2, and 
are scalars; w^^ and ^-re vectors; and is a (traceless) matrix. 

The excess free energy density in Eq. Q separates naturally into three parts: 

$ = $1 + $2 + ^3, (8) 
which are defined as 

"^1= E (9) 

i=C,P 

$2 = E (^>2 - nti • nQ (n^, nf ) , (10) 

iJ=C,P 

^3 = ^ E f^4«'2^2 - 4 ni2 ■ n^2 

[<2<2nv2 - tr (<2n42n^2)] ) </'ufc (^^ , <) • (H) 



+ 2 



Here tr denotes the trace operation and 

<^i...k{Vc,Vp) = ^pFo^ivcVp) (12) 

are derivatives of the zero-dimensional (Od) excess free energy, 

pFodirjc, Tjp) = {l-ric- rjp) ln(l - rjc) + rjc- (13) 
Equations (p?|)- |T!^ completely specify the excess free energy in Eq. (^). 

3.2. Application to Bulk Fluids 

For bulk fluid states, the density proflles are spatially constant. In the homogeneous 
limit [pc(r) — > pc and pp{t; Rp) — > pp{Rp)], Eq. ((H) yields the bulk fluid free energy 
density: 



V 



In (pcAc) PP^^P) {pp{Rp)Al) - 1 

+ <Phsivc) - / dRpppiRp) [\na{Rp; r]c) - /chain(i?p; Vc)] , (14) 
Jo 
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which is still a functional of the polymer density distribution, pp{Rp). Here 4>hs{vc) is 
the excess free energy density of the pure HS system, given by 

. . ^ 3Vc[Svc{2 - vc) - 2(1 - VcrH^ - Vc)] .... 

'^^^^''^^ SnRUl - VC? ' ^ ^ 

the same result as in the scaled-particle and Percus-Yevick compressibility 

approximations [301; and a{Rp;r]c) is the free volume fraction, i.e., the fraction of 
the total volume not excluded to the polymer by the HS colloids, given implicitly by 

In a{Rp- rjc) = ln(l - Vc) - C™7'", (16) 

m=l 

where 7 = Vc/{^ ~ Vc) and the coefficients are polynomials in the polymer-to-colloid 
size ratio, q = Rp/Rc- Ci = 3q + Sg^ + q^^ C2 = 9gV2 + 3q^, and C3 = 3q^. 

It is worth noting that the free volume fraction is related to the polymer one-particle 
direct correlation function, 

c?{r;Rp; [pc{r)]) = -B ^/"^ (17) 

dpp{r; Rp) 

which has its physical origin in colloid-polymer correlations. Substituting our DF 
approximation for F^^ into Eq. ()17|) and taking the homogeneous limit, we obtain 

^^\Rp- 7^c) = -J2q^* Wu{Rp) = In a(i?p; r/c), (18) 

where * denotes a convolution. Because our approximate excess free energy functional 
is linear in the polymer density, Cp\Rp;r]c) is independent of pp{Rp), depending only 
on the colloid density. 

In equilibrium, the polymer density distribution in Eq. (|T^ is fixed by the Euler- 
Lagrange equation for the polymers: 

|^^ = /ip(^p), (19) 
dpp{Rp) 

where fip{Rp) is the chemical potential of polymers with radius of gyration Rp. 
Substituting Eq. (fT^ into Eq. (fT^ . we have 

\n[pp{Rp;r]c)Al) - lnaiRp;r]c) + PU.URp) = PMRp). (20) 
where pp{Rp; rjc) is the equilibrium polymer density distribution at bulk colloid volume 
fraction rjc- Here we make the simplifying assumption that the chain free energy of 
a polymer with given Rp in the system is the same as that of an equal-sized polymer 
in the reservoir, and thus independent of colloid concentration. Because the system 
and reservoir must be in equilibrium with respect to polymer exchange, the right side 
of Eq. can be equated with the chemical potential of the corresponding polymer 
species in the reservoir: 

pip{Rp) = pi'p{Rp) = (3-^ In [p'-p{Rp)Kl) + U^^{Rp), (21) 

where p'^p{Rp) denotes the polymer density distribution in the reservoir. Further progress 
is facilitated by factoring the reservoir polymer density distribution, according to 

p^Rp) = p'pP^iRp), (22) 
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where pp and P^{Rp) are the mean density and radius-of-gyration probabihty 
distribution, respectively, of polymers in the reservoir. The advantage of Eq. (j^^ is 
that the probability distribution is simply related to the chain free energy via |2ZI 

/3 /chain (-Rp) = — lnP''(_Rp) + constant independent of Rp, (23) 

the arbitrary additive constant being chosen to normalize the distribution: 
j;^dRpP''{Rp) = 1. Combining Eqs. ^ and and using Eqs. ^ and 

we obtain 

ppiRp;Vc) = Ppa{Rp;r]c)P'{Rp) = a{Rp;vc)Pp{Rp). (24) 
The mean polymer densities in the system and reservoir are now seen to be related via 

ppiVc) = dRp pp{Rp;ric) = ppaesivc), (25) 
Jo 

where 

aes{vc)= dRpa{Rp;r]c)P''-{Rp) (26) 
Jo 

is an effective free volume fraction. Equation is a generalization of the 

incompressible-polymer relation, pp{rjc) = Ppa^rjc), from free volume theory [3 
Factoring pp{Rp]rjc) according to 

Pp{Rp: Vc) = Pp{^c)P{Rp; r]c), (27) 
serves to define 

P{Rp;Vc) = '^^^f^P'iRp) (28) 

as the normalized radius-of-gyration probability distribution of polymers in the system. 
Equation p8|l makes manifest that the radius-of-gyration distribution in the system 
depends on colloid concentration and differs from the distribution in the reservoir. 

Finally, substituting Eqs. (^Hj) and into Eq. (fT^ and rearranging, we obtain 
the equilibrium bulk fluid free energy density. 



In (pcAc) - l] + Pp^csivc) [HPp^D - l] + <PUvc). (29) 



Equations (j28|) and (|29|) are the main theoretical results of the paper and our basis 
for calculating polymer size distributions and phase behaviour in Sec. |3 Practical 
implementation, however, first requires specification of P^{Rp). 



3.3. Radius of Gyration Probability Distribution 

To approximate the radius-of-gyration probability distribution, we appeal to the 
statistical mechanics of a freely-jointed chain. The configuration of a chain of links, 
each link of length a, can be specified by a set of position vectors, {Ri} = (Ro • • • Rat), 
of the joints. For a given configuration, the radius of gyration, Rp, is defined by jSI] 

N 

Rp = ~ Rem) , (30) 

i=l 
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where Rem = Sili R-i is the centre-of-mass position of the chain. Averaged over 
configurations, the radius of gyration follows the probability distribution, P^{Rp). 
Moments of the distribution are defined as 

{R}>) = dRpR},P'{Rp), (31) 
Jo 

the second {n = 2) moment being related to the rms radius of gyration of polymers in 
the reservoir |SI1IS2], -Rg = \/ (Rp) = a\J N/Q. For polymers in the system, the rms 
radius of gyration is defined as 

1/2 



R^{t1c)= (yj^ dRpRlP{Rp-Tic) 



= r dRpRla{Rp-r^c)P'{Rp)] , (32) 

which clearly varies with bulk colloid concentration. Note that -Rg(?7c = 0) = Rg- 

In general, the radius of gyration is a more appropriate measure of polymer size than 
the end-to-end displacement, which is not well defined for branched polymers. However, 
in contrast to the simple Gaussian distribution of the end-to-end displacement, the 
radius-of-gyration distribution is nontrivial, even for ideal chains. Flory and Fisk [^ 1^ 
first proposed an empirical approximation for P^{Rp) based on the exact even moments 
calculated by Fixman jHSj. Subsequently, Fujita and Norisuye [321 EH] calculated the 
distribution exactly, obtaining the analytical result 



X I 1 - J-^ Ky^{tk) + fl - J-) ^3/4(4) 



?r\2 



(33) 



where t = (Rp/R^^) , are the modified Bessel functions of the second kind, and 
tk = (4A; + 3)^/(8t). To confirm Eq. we have carried out Monte Carlo simulations of 
an ideal chain for = 100 and 1000. Histograms of Rp, generated from 10^ independent 
configurations, are in essentially perfect agreement with the analytical expression |37j . 

Equation (j3^ specifies the radius-of-gyration probability distribution and thus, 
together with Eq. (j^^ . the chain free energy. We emphasize that in modelling the chains 
as freely jointed, and so applying the ideal-chain distribution [Eq. (j221)] to polymers 
mixed with colloids, we implicitly neglect any effect of colloidal confinement on the basic 
shape of the distribution. In particular, confinement-induced nonspherical distributions 
are neglected, although variation in the polymer size is allowed. 

4. Results: Demixing Phase Behaviour and Polymer Size Distribution 

The general conditions for phase coexistence are equality of the total pressure, 

(34) 

^ i=C,P ^1^1- 
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and of the chemical potentials for each species, 

^, = ^i|M, i = C,P, (35) 

between the coexisting phases. Equilibrium between phases I and II requires pl^^ = p"^ 
and /i^ = /i", equality of the polymer chemical potentials being enforced by Eq. (fT^ . 

In practice, bulk fluid phase diagrams can be computed as follows. The effective 
free volume fraction, aesivc), is first determined by substituting Eqs. (fTHjl and (jHSI) 
into Eq. (pUj). Then, for a given reservoir polymer density, pp, Eqs. (|^. (jMI), and 
fl35|l (z = C) are solved numerically for the coexisting colloid packing fractions, tjq 
and r]^. Finally, Eq. converts from reservoir to system representation, giving the 
corresponding system polymer packing fractions, r/p and r/p. 

Figures Q and El present bulk fluid phase diagrams for rms reservoir polymer-to- 
colloid size ratios EJg/Rc =0.5 and 1, in both the system and reservoir representations. 
For comparison, we include demixing binodals both from the present theory, which 
takes into account polymer compressibility, and from the free volume theory of the AO 
model with incompressible polymer. For sufficiently high colloid and polymer packing 
fractions, the system demixes into colloid-rich (liquid) and colloid-poor (vapour) phases. 
We disregard here the liquid-solid branch of the phase diagram, since our size ratios 
are sufficiently high that fiuid-fiuid demixing can be assumed stable. The effect of 
polymer compressibility evidently is to shift the colloidal liquid branch of the binodal 
toward higher polymer density, stabilizing the system against demixing. Interpreted 
in terms of effective depletion-induced attraction between pairs of colloids, polymer 
compression shortens the range of attraction, tending to favour mixing. In passing, it 
may be anticipated that additional degrees of freedom allowing for nonspherical polymer 
conformations would tend to further stabilize the mixture if, as an alternative to bulk 
demixing, the polymer may simply distort its shape. 

Figures El and 01 show the corresponding normalized polymer size distributions as a 
function of radius-of-gyration, scaled to the rms reservoir value. For given rms reservoir 
size ratio, with increasing colloid concentration, the distributions narrow and shift 
toward smaller radii, reflecting compression of the polymer by the confining colloids. 
In the case of Kg/Rc = 0.5, the polymer is compressed, relative to its rms size in the 
reservoir, to Rg/R^g =0.974, 0.938, 0.891, and 0.829 at colloid packing fractions rjc =0.1, 
0.2, 0.3, and 0.4, respectively. For R^^/Rc = 1, the polymer is even more strongly 
compressed to Rg/Rg =0.915, 0.826, 0.738, and 0.653 at the same respective values of 
rjc- Note that, for given rjc, P{Rp) does not depend on rjp (nor on rjp). Hence in the 
phase diagram, the polymer size distribution is constant along vertical lines (in both 
system and reservoir representations). What changes, of course, is the overall prefactor 
in pp{Rp). This invariance is an approximate feature of the theory resulting from the 
linearity of the excess free energy in polymer density. In reality, the presence of the 
polymers may alter the confining structure of the colloids, in turn modifying P{Rp). 
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5. Conclusions 

We have investigated the bulk fluid demixing behaviour of a model mixture of hard 
colloidal spheres and non- adsorbing polymer chains in a theta solvent. The polymer 
component is modelled, on a mesoscopic level, as a polydisperse mixture of effective 
spheres with radii free to adjust to allow for colloid- induced compression. The model 
goes beyond previous studies by treating the polymer radius of gyration as an internal 
degree of freedom that is polydisperse and varies between coexisting phases. Like all 
effective-sphere models, however, our approach ignores details of the polymers on the 
segment level and thereby neglects any effects of polymer shape anisotropy. 

To describe the model system, we have developed a geometry-based density- 
functional theory that incorporates polymer configurational free energy from 
consideration of the statistical mechanics of an ideal (random-walk) chain. For 
simplicity, the polymer chain free energy is assumed to be insensitive to colloidal 
confinement. Minimization of an approximate free energy functional with respect to 
the polymer size distribution yields both the equilibrium size distribution - modified 
by colloid-polymer interactions - and the free energy, from which we compute bulk 
fiuid phase diagrams. Polymer compression is found to increase with increasing colloid 
concentration and polymer-to-colloid size ratio. Correspondingly, the demixing fluid 
binodal shifts to higher polymer densities on the colloid-rich branch, favouring mixing. 

Predictions of the theory may be tested by experiment and simulation. Although 
the predicted shift of the demixing binodal due to polymer compression would appear 
to improve agreement with limited experimental data near the theta point 8j, more 
detailed measurements are certainly desirable. Polymer size distributions could be 
probed using either light scattering, by index matching the colloids (but not polymers) to 
the solvent, or neutron scattering by selectively deuterating the polymer (hydrocarbon) 
backbones. Detailed molecular simulations of segmented polymer chains mixed with 
hard spheres, along the lines of Refs. (23j and [21], could provide the cleanest tests of 
our approximations and predictions. 

Future applications to inhomogeneous colloid-polymer mixtures under the influence 
of external potentials (due to walls, gravity, etc.) and to liquid-solid transitions would 
be feasible and worthwhile. Extensions of the theory could explore the sensitivity 
of the polymer chain free energy to colloidal confinement, polydispersity in chain 
length, influences of anisotropic polymer conformations, distinctions between linear and 
branched polymers (as in colloid-star polymer mixtures jSHl IHH HO]), and effects of 
polymer nonideality in good and poor solvents. 
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Figure 1. Fluid- fluid demixing phase diagram for rms reservoir polymer-to-colloid size 
ratio Rg/Rc = 0.5. Solid curve: binodal from the present theory, taking into account 
polymer compressibility; dashed curve; binodal from free volume theory for the AO 
model with incompressible polymer. Thin straight lines indicate tie lines between 
coexisting phases; dots indicate critical points, (a) Reservoir representation: polymer 
reservoir packing fraction rjp vs. colloid packing fraction r]c; (b) system representation: 
system polymer packing fraction rjp vs. rjc- 
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Figure 3. Normalized probability distribution, P{Rp), of scaled polymer radius 
of gyration, Rp/R^ for rms reservoir polymer-to-colloid size ratio Rg/Rc = 0.5 
(corresponding to Fig. Colloid packing fraction, rjc, increases in direction of 
arrow: rjc =0, 0.1, 0.2, 0.3, and 0.4. Corresponding rms polymer size ratios decrease: 
Rg/Rl =1, 0.974, 0.938, 0.891, and 0.829. 
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Figure 4. Same as Fig. O but for R^/Rc ~ 1 (corresponding to Fig.j^J. In direction 
of arrow, rms polymer size ratios are Rg/Rg —I, 0.915, 0.826, 0.738, and 0.653. Note 
the change in scale compared with Fig. |31 



